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Abstract 



We consider one-dimensional chain of coupled linear and nonlinear oscillators 
with long-range power-wise interaction. The corresponding term in dynamical equa- 
tions is proportional to l/|n — m\ a+1 . It is shown that the equation of motion in 
the infrared limit can be transformed into the medium equation with the Riesz 
fractional derivative of order a, when < a < 2. We consider few models of 
coupled oscillators and show how their synchronization can appear as a result of 
^ ■ bifurcation, and how the corresponding solutions depend on a. The presence of 

fractional derivative leads also to the occurrence of localized structures. Particu- 
lar solutions for fractional time-dependent complex Ginzburg-Landau (or nonlinear 
Schrodinger) equation are derived. These solutions are interpreted as synchronized 
states and localized structures of the oscillatory medium. 
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Although the fractional calculus is known for more than two hundred years and 
its developing is an active area of mathematics, appearance and use of it in physical 
literature is fairly recent and sometimes is considered as exotic. In fact, there 
are many different areas where fractional equations, i.e., equations with fractional 
integro-differentiation, describe real processes. Between the most related areas are 
chaotic dynamics pQ, random walk in fractal space-time [2 J and random processes of 
Levy type |3| \Q \5[ |U] . The physical reasons for the appearance of fractional equations 
are intermittancy, dissipation, wave propagation in complex media, long memory, 
and others. This article deals with long-range interaction that can work in some 
way as a long memory. A unified approach to the origin of fractional dynamics 
from the long-range interaction of nonlinear oscillators or other objects permits 
to consider such phenomena as synchronization, breathers formation, space-time 
structures by the same formalism using new tools from the fractional calculus. 

1 Introduction 

Collective oscillation and synchronization are the fundamental phenomena in physics, chemistry, 
biology, and neuroscience, which are actively studied recently 0IB1IEJ) having both important 
theoretical and applied significance. Beginning with the pioneering contributions by Winfree ^01 
and Kuramoto studies of synchronization in populations of coupled oscillators becomes an 
active field of research in biology and chemistry. An oscillatory medium is an extended system, 
where each site (element) performs self-sustained oscillations. A good physical and chemical 
example is the oscillatory Belousov-Zhabotinsky reaction |121ll3llTT] in a medium where different 
sites can oscillate with different periods and phases. Typically, the reaction is accompanied by 
a color variation of the medium. Complex Ginzburg-Landau equation is canonical model for 
oscillatory systems with local coupling near Hopf bifurcation. Recently, Tanaka and Kuramoto 
|14j have shown how, in the vicinity of the bifurcation, the description of an array of nonlocally 
coupled oscillators can be reduced to the complex Ginzburg-Landau equation. In Ref. [T5] . 
a model of population of diffusively coupled oscillators with limit-cycles is described by the 
complex Ginzburg-Landau equation with nonlocal interaction. Nonlocal coupling is considered 
in Refs. ^01 El El- The long-range interaction that decreases as with < a < 2 
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is considered in Refs. |181 1191 I2U1 1211 I22| with respect to the system's thermodynamics and 
phase transition. It is also shown in [23] that using the Fourier transform and limit for the wave 
number k — > 0, the long-range term interaction leads under special conditions to the fractional 
dynamics. 

In the last decade it is found that many physical processes can be adequately described by 
equations that consist of derivatives of fractional order. In a fairly short period of time the list 
of such applications becomes long and the area of applications is broad. Even in a concise form, 
the applications include material sciences [22] ESI EH1 E2 , chaotic dynamics 1; , quantum theory 
[251 l%9l QUI l3*T] , physical kinetics ^ El E31 E2] > fluids and plasma physics jHU ES] , and many 
others physical topics related to wave propagation [35], long-range dissipation [37], anomalous 
diffusion and transport theory (see reviews ^ [2] [^J [3S] ) . 

It is known that the appearance of fractional derivatives in equations of motion can be linked 
to nonlocal properties of dynamics. Fractional Ginzburg-Landau equation has been suggested 
in |391 140[ l4"T] . In this paper, we consider the synchronization for oscillators with long-range 
interaction that in continuous limit leads to the fractional complex Ginzburg-Landau equation. 
We confirm the result obtained in J23 . that the infrared limit (wave number k — > 0) of an infinite 
chain of oscillators with the long-range interaction can be described by equations with fractional 
Riesz coordinate derivative of order a < 2. This result permits to apply different tools of the 
fractional calculus to the considered systems, and to interpret different system's features in an 
unified way. 

In Sec. 2, we consider a systems of oscillators with linear long-range interaction. For infrared 
behavior of the oscillatory medium, we obtain the equations that has coordinate derivatives of 
fractional order. In Sec. 3, some particular solutions are derived with a constant wave number 
for the fractional Ginzburg-Landau equation. These solutions are interpreted as synchroniza- 
tion in the oscillatory medium. In Sec. 4, we derive solutions of the fractional Ginzburg-Landau 
equation near a limit cycle. These solutions are interpreted as coherent structures in the os- 
cillatory medium with long-range interaction. In Sec. 5, we consider the nonlinear long-range 
interaction of oscillators and corresponding equations for the spin field. Finally, discussion of 
the results and conclusion are given in Sec. 6. In two appendices there are details of derivation 
of fractional phase (spin) equation. 
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2 Long-range interaction of oscillators 
2.1 Linear nonlocal coupling 

Let us consider the oscillators that are described by 

±z n (t) = F(z n ), (1) 

where z n is the position of the n-th oscillator in the complex plane, and F is a force. As an 
example, for the oscillators with a limit-cycle, F can be taken as 

F(z) = (1 + ia)z - (1 + ib)\z\ 2 z. (2) 

For a = b, each oscillator has a stable limit cycle at \z n \ = 1. Consider an infinite population 
(N — * oo) of oscillators ((Hi. with linear nonlocal coupling 

d 

dl' 



-z n (i) = F(z n ) + g J a (n - m)(z„ - z m ), (3) 



where the nonlocality is given by the power function 

J a (n) = \n\- a -\ (4) 
This coupling in the limit a — > oo is a nearest-neighbor interaction. 

2.2 Derivation of equation for continuous oscillatory medium 

In this section, we derive the equations for continuous medium that consists of oscillators (|3"|) 
with long-range interaction @. Let us define field 

i r+oo +oo 

Z(x,t) = — dk e ikx ^ e- ikn z n (t) (5) 

to describe the continuous oscillatory medium. 

Multiplying Eq. © on exp(— ikn), and summing over n from — oo to +oo, we obtain 

OO , CO oo 1 

S e " 4fcn ^ z " (t)= E ^ fcn ^)+^o E E e " %_ m .<H-i (*"-*»)• ( 6 ) 

n=— oo n=— oo n=— oo m^n 

It is convenient to introduce a new field variable. Let us define 

oo 

y(M)= £ e"^ n (t), (7) 



J a (k) = E e-^Un) = £ e- tkn T^TI- (8) 

The interaction term of © can be presented as two terms 

00 1 
^ e 4fcn , _ m i a+ i ( 2 « ~ 2 ™) = 



n — m 

n=-oora/n 



= E E e ~ lkn \ n - m \ a +i Zn - E Y. e ~ lkn \ n _ m \ a +i Zm - ( 9 ) 

n=—oom^n n=-oom/n 

For the first term in r.h.s. of @: 

00 1 00 1 

Er e -ifcn 1 z = e ~ikn z 1 = 

^ \n-m\ a+1 n ^ n ^\n-m\ a+1 

n=-oom^n n=— 00 m^n 



00 1 

= E e ' lkn ^ E u^n = rtM)£(o), (10) 

where 

00 

= E roi = 2 E u^TT = 2C(« + 1). (11) 

n^O 1 1 71=1 1 1 

and C( z ) is the Riemann zeta- function. For the second term in r.h.s. of @: 

1 ^ ^ 1 



E E e ifcn | n _ m i Q +i^- E z -E e „ - ,„i«+i 

n=—oom^n m=—co nj^m 

= Y, z me~ ikm E e ^OE+T = *(*.*)&(*). (I 2 ) 

where J Q (A;) is defined by ©• 
As the result, Eq. © is 

J^O, t) = F(y(k, t)) + g [J a (0) - J a (k)]y(k, t), (13) 

where F(y(k,t)) is an operator notation for the Fourier transform of F(z n ). The function J a (k) 
introduced in (JSJ) can be transformed as 

OO -, — OO 



n^O 71=1 n=— 1 



-ikn 





Figure 1: The function J a (k) for orders a = 1.1, and a = 1.9. 



E c -ifcn 1 _|_ c ikn 1 _ 1 



n=l 



n=l 



n=l 



g— ifcn _|_ gjfcn 



= Li a+1 (e ife )+LWi(e^ fc ), (14) 

where Li a (z) is a polylogarithm function. This presentation was also obtained in [231, and it 
plays an important role in the following transition to fractional dynamics of J a {k) through the 
polylogarithm. Using the expansion 



Li p (e z ) = ni-P){-zf^ + Y J - 



n=0 



IV. 



-z n , \z\ < 2tt, 



we obtain 

oo 

J a (k) = 2T(-a) cos(Tra/2)\k\ a + 2 ^ 
Using (|14|) in the form 



c{a+1 - 2n) i-kr, j a (o)=2 C (o + i) 



n=0 



J a {k) = 2^ 



(2n)! 
cos(/cn) 



n=l 



a+l 



we can see that 



J a (k + 2nm) = J a (k), 



(15) 



(16) 



(17) 



(18) 



where m is an integer. For a = 2, J a (fc) is the Clausen function Cl^(k) j^J. The plots of J a (k) 
for a = 1.1, and a = 1.9 are presented in Fig. 1. 
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Then Eq. (|13|) has the form 

^y(k, t) = F(y(k, t)) - goa a \k\ a y(k, t) - 2g C( " 2n) (-k 2 ) n y(k, t), (19) 

n=l ^ 

where 

a a = 2T(-a) cos(vra/2), (0 < a < 2, a / 1). (20) 
To derive the equation for field ®, we can use definition @ 



1 



Z(x,t) = — e tkx y(k,t)dk, (21) 
and the connection between Riesz fractional derivative and its Fourier transform |42| : 

fc2 "-^' (22> 

The properties of the Riesz derivative can be found in |421 1431 El I45j . Its another expression is 
(> '" Z ^ f ) = ~ o^L 7o\ inZ(x,t) + V a _Z(x,t)) , (23) 



d\x\ a v ' ' 2cos(?ra/2) 
where a ^ 1,3,5..., and T>± are Riemann-Liouville left and right fractional derivatives 

1 d n f x Z(£,t)d£ 



1 d n f 

V%Z(x,t) = - / 

1 (n — a) ox J_ 



a— n+ 1 ' 

— oo ■ ' 



(x - C 



Substitution of Eqs. into Eq. (|2*3*|) gives 

9|x|« V ' J 2cos(vra/2)r(n-a)9x" Vi-oo (^-e) Q " n+1 A - x) a ~ n+1 J ' { ' 
Multiplying Eq. (|19() on exp(ifcc), and integrating over k from — oo to +oo, we obtain 

71=1 

The first term (n = 1) of the sum is ((a — l)d 2 Z. Let us compare the coefficients of terms 
with fractional and second derivatives in Eq. 1)26(1. For a — * 2, one can use the asymptotics 



C(«-l)~— ^ + 0(1), a a ^^— + 0(l) (a/2) 
a — 2 a — 2 

As an example, for a = 1.99, 

C(a - 1) » -99.42351, a Q « -100.92921. 
Therefore C(a - l)/a a ~ 1 for 2 - a < 1. 



2.3 Infrared approximation 

In this section, we derive the main relation that permits to transfer the system of discrete 
oscillators into fractional differential equation. This transform will be called infrared limit. For 
< a < 2, a/1 and k — > 0, the fractional power of \k\ is a leading asymptotic term in Eq. 
Ifigj). and 

[J a (0) - J a (k)] Rs a a \k\ a , (0<a<2, a^l). (27) 

Eq. (|27|) can be considered as an infrared approximation of (|19|) . Substitution of (|27|) into 
gives 

J^y(M) = %(M)) - Sba«l*f tf(M), (0<a<2, a^l). (28) 

Then 

^-Z = F(Z)+g a a -^Z, (0 < a < 2, a ^ 1), (29) 

where Z = Z(x,t) is defined by (|21|) . Eq. ()29j) can be considered as an equation for continuous 
oscillatory medium with a < 2 in infrared (k — > 0) approximation. 

As an example, for F(z) = Eq. Q29[) gives the fractional kinetic equation: 

^-Z = go a a ^-^Z, (0<a<2, a ^ 1) (30) 

that describes the fractional super diffusion 001115. For F{z) defined by ©, Eq. (|2*§)) is a 
fractional Ginzburg-Landau equation that has been suggested in [3^1 (see also |H3I^), and will 
be considered in Section 3. For a > 2 and k — > the main term in ()16|) is proportional to k 2 and 
in (|29|). H30|) . we have second derivative instead of the fractional one. The existence of critical 
value a = 2 was obtained in [23]. 

3 Fractional Ginzburg-Landau (FGL) equation 

3.1 Synchronized states for Ginzburg-Landau equation 

The one-dimensional lattice of weakly coupled nonlinear oscillators is described by 

^n(t) = (1 + ia)z n - (1 + ib)\z n \ 2 z n + (ci + ic 2 )(zn+i - 2z n + z n -i), (31) 

where we assume that all oscillators have the same parameters. A transition to the continuous 
medium assumes that the difference z n+ \ — z n is of order Ax, and the interaction constants 



ci and C2 are large. Setting c\ = g(Ax) 2 , C2 = gc(Ax) 2 , and Z(x,t) « Z(nAx,t) = z n (t), we 
get 

— Z = (1 + ia)Z - (1 + ib)\Z\ 2 Z + 5 (1 + ic)-^Z, (32) 

which is a complex time-dependent Ginzburg-Landau equation. The simplest coherent structures 
for this equation are plane- wave solutions 8 : 

Z(x,t) = R{K)exp[iKx - iu(K)t + 0„] , (33) 

where 

R{K) = {l-gK 2 ) l l\ u;(K) = (b-a) + (c-b)gK 2 , (34) 
and #o is an arbitrary constant phase. These solutions exist for 

gK 2 < 1. (35) 

Solution (|3*3*|) can be interpreted as a synchronized state [SJ- 

3.2 Particular solution for FGL equation 

Let us come back to the equation for nonlinear oscillators (J2J) with F(z) in Eq. © and long-range 
coupling (@J: 

^-z n = (1 + ia)z n - (1 + i6)|z n | 2 z„ + 5o ^ i — _ 1 On ~ z m)> (36) 

mjtn 

where z n = z n (i) is the position of the n-th oscillator in the complex plane, 1 < a < 2. The 
corresponding equation in the continuous limit and infrared approximation can be obtained in 
the same way as (j29jl : 



^-Z = (1 + io)Z - (1 + ib)\Z\ 2 Z + 5 (1 + ic)J^Z, (37) 



where (7(1 + ic) = goa a , and 1 < a < 2. Eq. (|37jl is a fractional generalization of complex time- 
dependent Ginzburg-Landau equation (|3*2l (compare to ()29[)). Here this equation is derived in 
a specific approximation for the oscillatory medium. 
We seek a particular solution of (|37|) in the form 



Z(x,t) = A(K,t)e iKx , (38) 
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which allows to use 

iKx I j^\a iKx 



K\ a e lKx . (39) 



d\x\ a 

Eq. ("""""[) represents a particular solution of ((57)1 with a fixed wave number if. 
Substitution of ((3*H)l into ((3*7)) gives 

^A(K, t) = (l + ia)A - (1 + ib)\A\ 2 A - g{l + ic)|if | a A (40) 

Rewriting this equation in polar coordinates, 

A{K,t) = R{K,t)e^ K *\ (41) 

we obtain 

^ = {l-g\K\«)R-&, 

^- = {a-cg\K\ a )-bR 2 . (42) 
The limit cycle here is a circle with the radius 

R = (1 -£|i\T) 1/2 , £/|ET| a < 1. (43) 

Solution of 1(42)1 with arbitrary initial conditions 

i?(if,0) = i? , 0(tf,O)=0 o (44) 



is 



R(t) = Roil - g \K\"?l* (Rl + {l-g\K\«- ^) e -ad-^')*) 1/2 , (45) 



0(t) =-|ln[(l- IT 1 («o + (1 " <?I*T " ^o)e- 2a *)] " "a(K)t + O> (46) 

where 

u a (K) = (b-a) + (c-b)g\K\ a , 1 - g\K\ a > 0. (47) 

This solution can be interpreted as a coherent structure in nonlinear oscillatory medium with 
long-range interaction. 



If 

?2 



R z = l-g\K\ a , g\K\ a <l, 
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then Eqs. (|4"5j) and (HHJ) give 

R(t) = Ro, 6(t) = -u a (K)t + 9 . (48) 

Solution Q48[) means that on the limit cycle (|43|) the angle variable rotates with a constant 
velocity a; a (iC). As the result, we have the plane- wave solution 

Z(x, t) = (l- g |if|«)i/2 e ^-iu, a (X)t+i0 o 5 i _ > o, (49) 

which can be interpreted as synchronized state of the oscillatory medium. 

For initial amplitude that deviates from ()43j) . i.e., / 1 - gjiTl , an additional phase 
shift occurs due to the term which is proportional to b in ()46|) . The oscillatory medium can be 
characterized by a single generalized phase variable. To define it, let us rewrite (|42|) as 

±lnR=(l-g\Kn-R 2 , (50) 
at 

j t = (a-cg\K\ a ) -bR 2 . (51) 
Substitution of R 2 from (|5U)) into (|5*T|) gives 

j t (9-b\nR) = (a-cg\K\ a )-b{l-g\K\ a ). (52) 

Thus, the generalized phase jS] can be defined by 

<p{R,9) = 9 -blnR. (53) 

From (|52j) . we get 

±-cj> = -u a {K). (54) 

This equation means that generalized phase </>(-R, 0) rotates uniformly with constant velocity. 
For g|-ftT| a = (b — a)/(b — c) < 1, we have the lines of constant generalized phase. On (R,9) 
plane these lines are logarithmic spirals 9 — bhiR = const. The decrease of a corresponds to 
the increase of K. For the case 6 = instead of spirals we have straight lines <f> = 9. 

3.3 Group and phase velocity of plane-waves 



Energy propagation can be characterized by the group velocity 

dio a (K) 



Ja > 9 ~ dK 
11 



(55) 



From Eq. (|47|). we obtain 

v a , g = a(c-b)g\K\ a - 1 . (56) 



For 



we get 



The phase velocity is 



For 



\K\ < K x = (a/2) 2 ~ a , (57) 

\v a ,g\ > \v2,g\- (58) 

v atPh = u a (K)/K = (c-b)g\K\ a -\ (59) 

K <K 2 = 2 a ~ 2 , (60) 



we have 



ba,ph| > \v2,ph\- (61) 

Therefore long-range interaction decreasing as with 1 < a < 2 leads to increase the 

group and phase velocities for small wave numbers (K — > 0). Note that the ratio v a ^g/v a ^pji 
between the group and phase velocities of plane waves is equal to a. 

3.4 Stability of plane-wave solution 

Solution of (fl9|) can be presented as 

X = R(K,t)cos(9(K,t) + Kx), Y = R(K,t) sm(6{K,t) + Kx), (62) 

where X = X(K,t) = ReZ(x,t) and Y = Y(K,t) = ImZ(x,t), and R(K,t) and 6{K,t) are 
defined by ((H)) and (|46|). For the plane- waves 

X (x, t) = (l- g\K\ a ) 1/2 cos {Kx - uj a {K)t + 0„) , 



F (a;,i) = (1 -5|i^r) 1/2 sin(iTx- uj a {K)t + 6 Q ) , l- 5 |ir| a >0. (63) 

Not all of the plane-waves are stable. To obtain the stability condition, consider the variation 
of (HOI) near solution (|63|) : 

j t $X{K, t) = A U 5X + A 12 5Y, j^Y{K, t) = A 21 5X + A 22 5Y, (64) 

12 



where 5X and 5Y are small variations of X and Y, and 

A u = 1 - <?|i^r - 2X (X - bY ) - (X 2 + Y 2 ), 

A n = -a + gc\K\ a - 2Y (X - 6Y ) + 6(X 2 + Y 2 ), 
4u = a - gc\K\ a - 2X (Y + bX ) - 6(X 2 + Y 2 ), 

A 22 = l- g\K\ a - 2Y (Y + bX ) - (X 2 + Y 2 ). (65) 
The conditions of asymptotic stability for (|64[) are 

An + A 22 < 0, A11A22 - A 12 A 21 < 0. (66) 

From Eqs. (|63|) and (|65j) . we get 

A n + A 22 = -2(1 - g\K\ a ), l-g\K\ a >0, (67) 

and the first condition of (|6U)) is valid. Substitution of Eqs. and (|B31) into Q66|) gives 

A\\A 22 - A 12 A 21 = (6(1 - - (a - gc\K\ a )) (36(1 - 5 |ET| a ) - (a - gc\K\ a )) . (68) 

Then the second condition of ()66|) has the form 

(V- l)(V-3) < 0, (69) 

where 

As the result, we obtain 

< 1 - g\K\ a < a/b - {c/b)g\K\ a < 3(1 - g\K\ a ), (71) 

i.e. the plane- wave solution (|49|) is stable if parameters a, 6, c and g satisfy (|71jl. Condition ()71|) 
defines the region of parameters for plane waves where the synchronization exists. 
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3.5 Forced FGL equation for isochronous case 

In this section, we consider FGL equation (|4U|) forced by a constant E (the so-called forced 
isochronous case (b = 0) 8 ): 

— A = (l + ia)A-\A\ 2 A-g(l + ic)\K\ a A-iE, (ImE = 0) (72) 
at 

where A = A(K,t), and we put for simplicity 6 = 0, and K is a fixed wave number. Our 
main goal will be transition to a synchronized states and its dependence on the order a of the 
long-range interaction. The system of real equations is 

—X = (1 - g\K\ a )X - (a - gc\K\ a )Y - (X 2 + Y 2 )X, 
dt 

jY = (1 - g\K\«)Y + (a- gc\K\ a )X - (X 2 + Y 2 )Y - E, (73) 

where X = X(K,t) is real and Y = Y(K,t) imaginary parts of A(K,t). 

Numerical solution of Eq. (|73|) was performed with parameters a = 1, g = 1, c = 70, E = 0.9, 
K = 0.1, for a within interval a G (1; 2). The results are presented on Fig. 2, and Fig. 3. For 
cto < a < 2, where ao ~ 1.51..., the only stable solution is a stable fixed point. This is region 
of perfect synchronization (phase locking), where the synchronous oscillations have a constant 
amplitude and a constant phase shift with respect to the external force. For a < ao the global 
attractor for ()73() is a limit cycle. Here the motion of the forced system is quasiperiodic. For 
a = 2 there is a stable node. When a decreases, the stable mode transfers into a stable focus. 
At the transition point it loses stability, and a stable limit cycle appears. As the result, we have 
that the decrease of order a from 2 to 1 leads to the loss of synchronization (see Figures 2-3). 

In Fig. 2, (a = 2.00, a = 1.70, and a = 1.60, a = 1.56), we see that in the synchronization 
region all trajectories are attracted to a stable node. 

In Fig. 3, (a = 1.54, a = 1.52, and a = 1.50, a = 1.40), a stable limit cycle appears via 
the Hopf bifurcation. For a = 1.54, and a = 1.52, near the boundary of synchronization the 
fixed point is a focus. For a = 1.4, the amplitude of the limit cycle grows, and synchronization 
breaks down. 
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Figure 2: Approaching to the bifurcation point a = «o = 1.51... of solution of forced 
FGL equation for isochronous case with fixed wave number K = 0.1 is represented by real 
X(K, t) and imaginary Y(K, t) parts of A(K, t). The plots for orders a = 2.00, a = 1.70, 
a = 1.60, a = 1.56. 
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Figure 3: Transformation to the limit cycle of solution of forced FGL equation for 
isochronous case with fixed wave number K = 0.1 is represented by real X(K,t) and 
imaginary Y(K,t) parts of A(K,t). The plots for orders a = 1.54, a = 1.52, a = 1.50, 
a = 1.40. 
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3.6 Phase and amplitude for forced FGL equation 

The oscillator medium can be characterized by a single generalized phase variable (|53|) . We can 
rewrite (|53")l as 

Y) = arctan(Y/X) - - ln(X 2 + Y 2 ), (74) 
where X and Y are defined by (j62j) . For E = 0, the phase rotates uniformly 

^-6 = -u> a {K) = a-gc\K\ a , (75) 
at 

where cj a (if) is gived by (|47|) with 6 = 0, and can be considered as a frequency of natural 
oscillations. For E ^ 0, Eqs. (JTHJ) and (J71J) give 

^0 = -w a (ir) -Ecoscf). (76) 

This equation has an integral of motion. The integral is 

h = 2(cj 2 - J B 2 )- 1/2 arctan((cj-S)(cu 2 -^ 2 )- 1/2 tan(^(t)/2)) + t, J 2 > E 2 , (77) 

/ 2 = 2(S 2 -cj 2 )~ 1 / 2 arctanh^-u;)(S 2 -u; 2 )- 1 / 2 tan((/)(t)/2)) + 1, J 2 < E 2 . (78) 

These expressions help to obtain the solution in form ()41|) for forced case ()72(l keeping the same 
notations as in (|41j) . For polar coordinates we get 

^ = (1 - g\K\ a )R - R 3 -Esm9, 
dt y y| 1 ; 

^ = (a-c 5 |K| (79) 

Numerical solution of (|79|) was performed with the same parameters as for Eq. (|7cS|) . i.e., a = 1, 
g = 1, c= 70, £7 = 0.9, -fT = 0.1, and a within interval a £ (1, 2). The results are presented in 
Fig. 4 and Fig. 5. 

The time evolition of phase B(K,t) is given in Fig. 4 for a = 2.00, a = 1.50, a = 1.47, 
a = 1.44, a = 1.40, a = 1.30, a = 1.20, a = 1.10. The decrease of a from 2 to 1 leads to the 
oscillations of the phase 9(K,t) after the Hopf bifurcation at ao = 1.51..., then the amplitude 
of phase oscillation decreases and the velocity of phase rotations increases. 

The amplitude R(K,t) is shown on Fig. 5 for a = 1.6, a = 1.55, a = 1.55, a = 1.51, 
a = 1.50, a = 1.45, a = 1.2. The appearance of oscillations in the plots means the loss of 
synchronization. 
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Figure 4: Phase 6(K,t) for K = 0.1 and a = 2.00, a = 1.50, a = 1.47, a = 1.44, 
a = 1.40, a = 1.30, a = 1.20, a = 1.10. The decrease of order a corresponds to the 
clockwise rotation of curves. For upper curve a = 2. For the most vertical curve a = 1.1. 
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Figure 5: Amplitude R(K,t). The upper curve corresponds to a = 2 for all plots. The 
lower curves correspond to a = 1.6, a = 1.55, a = 1.51, a = 1.50, a = 1.45, a = 1.2. The 
appearance of oscillations on the plots means the loss of synchronization. 
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4 Space-structures from FGL equation 

In previous sections, we considered mainly time-evolution and "time-structures" as solutions for 
the FGL equation. Particularly, synchronization process was an example of the solution that 
converged to a time-coherent structure. Here we focuse on the space structures for the solution 
of FGL equation (|37|) with b = c = and the constants a\ and a 2 ahead of linear term: 

^-Z = {a l +ia 2 )Z-\Z\ 2 Z + g^—Z. (80) 
at o\x\ a 

Let us seek a particular solution of (|80jl in the form 

Z(x,t) = R(x,t)e w{t) , R*(x,t) = R(x,t), 9*(t) = 6{t). (81) 
Substitution of ((ED into (JHOJ) gives 

—R + iR-Q t m = («i + io, 2 )R -R 3 + 9q^R, (82) 



or 

§i R = "i R - R3 -^ R ' iw-* < 83 > 

1/2 

Using 6(t) = a 2 t + 6(0), we arrive to the existence of a limit cycle with Rq = a x . 
Let us find a particular solution of (|83|) in the vicinity of the limit cycle, i.e. 

R(x,t) = R + eR 1 +e 2 R 2 + (e < 1). (84) 

1 /2 

Zero approximation i?o = a/ satisfies (j83(l since 

1 = 0, (85) 



d\x 

and for R\ = i?i(x,t), we have 

d d a 

—Ri = -2a 1 R 1 + Oj^Ri ■ 86 
at o\x\ a 

Consider the Cauchy problem for (|8b[) with initial condition 

i2i(x,0) = ¥j(x). (87) 

By solution of the Cauchy problem we mean a function R±(x, t) which satisfies ()86|) and condition 
(|87|). The Green function (or fundamental solution) of the Cauchy problem is the (generalized) 
function G(x, t) such that 

/+oo 
G(x',t)(p(x - x')dx' 
-oo 
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The Green function is the solution of (|86|) with (p(x) = 5(x). 

Let us apply Laplace transform for t and Fourier transform for x: 

"OO f + OO 

-st+ikx/ 



G(k,s)= dt dx e- st+tkx G(x,t). (89) 

Jo J-oo 



By the definition of Riesz derivative, 

& 



-G(x,t)^-\k\«G(k,s), (90) 
and for the Laplace transform with respect to time 

^G(x, t) < — ► sG(k, s) - 1. (91) 
Applying (|55 fl - <fgT |> to (fSHjl. we obtain 

sG(k, s)-l = -2aiG(k, s) - g\k\ a G(k, s), (92) 



or 



G(k,s) = fTT-. (93) 

1 ' s + 2ai + g\k\ a V ; 



Let us first invert the Laplace transform in ()93|) : 

1 



< — ► e~ at , (94) 



s + a 

Then, the Fourier transform of the Green function 

/+oo 
dxe ikx G(x,t) (95) 
-oo 

has the form 

G(k, t) = e -(2ai+flW a )t = e -2a 1 t e - fl [fc|«t_ (96) 

As the result, we get 

G(x,t) = (gtr^e-^L^xigt)- 1 ^). (97) 



where 



i r+oo 

L a (x) = — dk e -^x e -a\k\" (gg) 

2tt 7_oo 



is a Levy stable p.d.f. [22] ■ The p.d.f. L a (x) for a = 2.0, a = 1.6, and a = 1.0 are shown on 
Fie. 6. 
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As an example, for a = 1 we have the Cauchy distribution with respect to coordinate 

e-^^L 1 (x) = -^—, (99) 

7T X z + 1 



and 



1 (at)~ l e~ 2ait 

G M = ;^iPTr (100) 



For a = 2, we get the Gauss distribution: 

T, n (rr\ — _ 

and 



e" fc2 «— ► L2(x) = -^e" 1 ' 2 / 4 , (101) 



G(x,i) = (jiJ-V^-^'i-e^/M. (102) 

2vvr 

For 1 < a < 2 the function L a (x) can be presented as the convergent expansion 

LJx) = ?(-x) n K , ' ' sin nvr/2). (103 

irx ^— ' n! 

n=l 

The asymptotic (x — ► oo, 1 < a < 2) is given by 

MaO y(-l) n x- n " r ^ 1+ , nQ) sin(mr/2), x -» oo, (104) 

7tx ^— ' n\ 

n=l 

with the leading term: 

L a (x) ~vr" 1 r(l + a)x" a ~ 1 , x^oo. (105) 
As the result, the solution of (|%U|) is 

Z(x, t) = e ^*+ e (°)) ^ a } /2 + e(^)-V« e - 2a i* L a {x> \gt)- l ' a ) tp{x - x')dx' + 0{e 2 )^j . 

(106) 

This solution can be considered as a space-time synchronization in the oscillatory medium with 
long-range interaction decreasing as |x| _( - a+1 - ) . 

For (p(x) = 8{x — xo), solution ()1U6|) has the form 

Z( x , t) = e i(a2t+m) [a] 12 + e{gt)- l l a e- 2ait L a {(x - x )(gt)- 1 / a ) + 0(e 2 )) , (107) 

and the asymptotic is 

Z(x, t) = e i( - a2t+e ^ (a{ /2 + e 5 te- 2ai *7r- 1 r(l + a)(x - xo)""" 1 + 0(e 2 )) , x -» oo. (108) 
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This solution shows that the long-wave modes approach to the limit cycle exponentially with 
time. For t = l/(2a\), we have the maximum of with respect to time: 

max\Z{x,t)\ = a/ 2 + eg^-^(x - x )- a ^ + 0{s 2 ). (109) 
*>o Z7re 

As the result, we have the power law decay with respect to coordinate for the space structures 
near the limit cycle \Z\ = a^ 2 . 



5 Nonlinear long-range interaction and fractional phase 
equation 

In this section, we would like to show one more application of the replacement of dynamical 
equation by the fractional ones for a chain with long-range interaction. The model was first 
considered in jlOl 111! I46j with application in biology and chemistry. This model has additional 
interest since it can be reduced to a chain of interacting spins. In this section, we will derive 
the fractional phase (spin) equations and leave their solution for the future publication. 

5.1 Nonlinear nonlocal phase coupling 

Let us consider the phase equation 

j£0n(t) =uJ n + 9^2 J a (n - m) sin(6 n - 9 m ), (110) 

where 6 n denotes the phase of the n-th oscillator, uj n its natural frequency, and 

J a (n) = In]- 01 ' 1 . (Ill) 

For a = — 1, Eq. ()110|) defines the Kuramoto model |llll46lH7| with sinusoidal non-local coupling 
(infinite radius of interaction) . We can rewrite Eq. (jllOj) for classical spin- like variables 

Sn (i) = e ^W. (112) 

Using 

s" 1 = s* n , 9 n = -ilns n , (113) 
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Figure 6: Gauss p.d.f. (a = 2), Levy p.d.f. (a = 1.6), and Cauchy p.d.f. (a = 1.0). Levy 
for a = 1.6 lies between Cauchy and Gauss p.d.f. In the asymptotic x — > oo and x > 3 
on the plot, the upper curve is Cauchy p.d.f, the lower curve is Gauss p.d.f. 
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we get 

j t 6 n = -ij t \ns n = -*s- l j t s n = -<j t *n, (H4) 

and 

sin(0 n -9 m ) = ^ ( e *(ffn-fl«) + e -Ko n -e m )\ = 1 ^ + <Sm) (n5) 
2z V / 2,1 



Then Eq. (fTTUl) is 



,— a n = iu> n + I E - +1 [s n s^ + S*g m ]. (116) 



This equation describes the long-range interaction of spin variables. We also will call Eq. 1)1161) as 
the phase coupling equation since \s n \ — const. Thermodynamics of the model of classical spins 
with long-range interactions (see for discussion Appendix 1) has been studied more than thirty 
years. An infinite one-dimensional Ising model with long-range interactions was considered 
by Dyson ^H]- The (i-dimensional classical Heisenberg model with long-range interaction is 
described in Refs. |19[ I20|. and its quantum generalization with long-range interaction decreasing 
as \n\~ a can be found in |21| . 

5.2 Phase-coupled oscillatory medium with nonlinear long-range 
interaction 

Let us derive an equation for continuous medium that consists of oscillators of (|11(J|) or H116|) 
type with nonlinear long-range interaction. The medium can be defined by the field 

S(x,t) = — dke ikx J2 e- ifcn s n (i)- (117) 

J 00 n=-oo 

We also will need the following momentum representations 

oo 

a(M)= E e " ,hs «W' ( 118 ) 

n=— oo 

For left hand side of (|116[) . we get (see Appendix 2): 

1 f°° 00 d H 

— dke lkx E e- lkn S * n -s n = S*(x,t)-S(x,t). (119) 

n=-oo 

For the interaction term, we obtain (Appendix 2): 
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(r) a °° tin -A- 1 — 2n) r) 2n \ 

2((a + l)S(x, t) - a a ^S(x, t) + 2 E {2n)[ Q^ S ^ *) ) > ( 120 ) 

where we use (fTo]) for Jq,(A;) which is defined by ©, and 

a a = 2r(-a) cos(vra/2). (121) 

For the term uj n , we use 



1 r+00 00 

- / dke ikx e~ ikn u n - (122) 



V ' 2vr 

If all oscillators have the same natural frequency u n = to, then u(x) = u>. 
As the result, 

d 

S*(x,t)—S{x,t) = iuj(x) - f a S*(x,t)S(x,t)- 



~9a (s*(x,t)^S(x,t) + S(x,t)-^S*(x,t)^ + 
+9 £ 1 | 2n)! J (** (*> 0^^, + S(x, t) ^S* (x, t) J , (123) 



where 



/« = 2 5 C(a + 1), = (l/2)a aS = ^r(-a) cos(to/2). (124) 

Eq. ()123j) is a fractional equation for oscillatory medium with long-range interacting spins (|116|) . 
We can call (|123|) the fractional phase equation. 

In the infrared approximation {k — > 0), we can use (|16jl : 

J a (fc) Rs 2r(-a)cos(vra/2)|/c| a + 2C(a + 1), < a < 2, a ^ 1, (125) 



and equation (|123j) is reduced to 

S*(x,t)^- t S(x,t)=iu(x)-f a -g a (s*(x,t)^S(x,t) + S(x,t)-^S*(x,t)), (126) 
where < a < 2, a ^ 1. 
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6 Conclusion 



One-dimensional chain of interacting objects, say oscillators, can be considered as a benchmark 
for numerous applications in physics, chemistry, biology, etc. All considered models were re- 
lated mainly to the oscillating objects with long-range power- wise interaction, i.e., with forces 
proportional to l/\n — m\ s and 2 < s < 3. A remarkable feature of this interaction is the 
existence of a transform operator that replaces the set of coupled individual oscillator equations 
into the continuous medium equation with fractional space derivative of order a = s — 1, where 
< a < 2, a ^ 1. Such transformation is an approximation and it appears in the infrared limit 
for wave number k —* 0. This limit helps to consider different models and related phenomena 
in a unified way applying different tools of fractional calculus. 

A nontrivial example of general property of fractional linear equation is its solution with a 
power-wise decay along the space coordinate. From the physical point of view that means a 
new type of space structures or coherent structures. The scheme of equations with fractional 
derivatives includes either effect of synchronization [Sj , breathers |531 154j , fractional kinetics , 
and others. 

Discrete breathers are periodic space-localized oscillations that arise in discrete and con- 
tinuous nonlinear systems. Their existence was proven in Ref. jSS]. Discrete breathers have 
been widely studied in systems with short-range interactions (for a review, see |571 [53 ). En- 
ergy and decay properties of discrete breathers in systems with long-range interactions have 
also been studied in the framework of the Klein-Gordon |54[ I58j . and the discrete nonlinear 
Schrodinger equations [5Ji]. Therefore, it is interesting to consider breathers solution in systems 
with long-range interactions in infrared approximation. 

We also assume that the obtained transform operator can be used for improvement of dif- 
ferent scheme of simulations for equations with fractional derivatives. 
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Appendix 1: Spin chain models with long-range inter- 
actions 

The model of classical spins with long-range interactions has been studied for more than thirty 
years. An infinite one-dimensional Ising model with long-range interactions was considered by 
Dyson ^B]- Existence of a phase-transition was proved for an infinite linear chain of spins s n 
with an interaction Hamiltonian 

H int = - J a (n-m)s n s m , J a {n) = \n\~ a ~ 1 . (127) 

It was proved that the Ising model (|127|) has a phase transition of the first order when < a < 1 , 
i.e. for fairly distant interactions. 

For the d-dimensional {d = 1,2) classical Heisenberg model with long-range interaction 

Hint = ~ ^2 J a{n-m) s n s m , (128) 

n, m^n 

the existence of the first order phase transition for d < a < 2d and nonexistence for a > 2d was 
proved in Refs. [HI EH- 

The quantum XY-model with long-range interaction (|128|) was considered in Refs. [2^. This 
model in the limit a — > oo corresponds to nearest-neighbor interaction. The spin operators s n 
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at each lattice site can be replaced by boson creation and annihilation operators: 

s* + isl = V2a+, s x n -is y n = V2a n . 

The site representation can be replaced by the momentum representation using Fourier trans- 
formation 

v N k ^ N k 

where a + (k), a{k) are the corresponding Fourier transforms of operators a+, a n , and N is a 
number of particles. Then the interaction part of the Hamiltonian can be rewritten as 

H int = Y J E{k)a + {k)a{k) (129) 

k 

with the dispersion relation 

N 

E{k) = [J o (0) - J a (k)], J a (k) = 2 J<x{n) cos(fcn). (130) 

71=1 

In the thermodynamic limit N — > oo the dispersion relation (jl.SOj) becomes 

n=l 

For small k (infrared limit \k\ — > 0) and < a < 2, a ^ 1, we have 

S(fc)«ja|fcr, (132) 



where 



7T9 



r(a) cos(7ra/2) 

For oscillatory medium, Eq. (|132|) gives the fractional Riesz derivatives of order a (compare 
to (HU». 

Appendix 2: Derivation of fractional phase equation 

The left hand side of (FTTH1) is 



dt " (2tt) 



1 /-+ 00 j 

/ dk 1 dk 2 e- mk 'e ink2 a*{k 1 ,t)-^a{k 2 ,t), (134) 

7>T J-oo «* 
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or for its Fourier transform 



e ~ lkns nJ t s n = T^yl / dhdk 2 e m ^- k ^a*(k 1 ,t)-a(k 2 ,t)= (135) 

i=— co ' ^ ' 00 n=— 00 

1 f°° d 
= — dk 1 dk 2 5(k 2 - h - k)a*(k 1 ,t)—a(k 2 ,t) = 
2tt J_ oq dt 

If 00 d 
= —J dk 2 a*{k 2 -k,t)-a{k 2 ,t). (136) 



Then 



1 f°° °° d 1 f'°° 1 f°° d 

— / dke ikx Yl e ~' tkns nT + s ™ = V / dke ihx — / dk 2 a*(k 2 - k,t)-a(k 2 ,t) = 

1 f 00 r/ 1 f'°° 

= ^J dk 2J t <k 2 ,t)—J dka*(k 2 -k,t)e ikx = 

1 f'°° d 1 f 00 d 

= — dk 2 e^ x -a(k 2 ,t) — dkm'farfe-* 1 * = S*(x,t)-S(x,t), (137) 
2ir J_ oc 27T J_ 00 dt 



-00 </ — 00 

where we use 

1 



+00 



S( x ,t) = — e lkx a(k,t)dk. (138) 
27T J_ 00 



Let us consider the interaction terms 



1 1 1 roo 

E |n-mh+i flm( * )fl " (t) = E (W I dk 1 dk 2 ^e^a{k,, t)a*(k 2 , t). 

(139) 

Multiplying this equation by exp(— ikn), and summing over n from — oo to +00, we obtain 



00 



E e -ikn \ 

n=—oo mj^n 



,«+! S n s m 

n — m \ a+1 



"~ 11 poo 

E ^ E | n _ m |« + i (27T) 2 / dk ld k 2 e^e- k Mh,t)a*(k 2 ,t) = 
1 r°° °° i 

= (2^/ ^^omkom) E E i re _ m |. + i etmfcle " m(fc+fc2) = 

^ ' 00 n=-oo m ^n 1 1 



I poo 00 ( 

/ dfcxdM*!, *)«•(*>,*) E ^ m(fc+fc2) E | n - m |a + l « <m/ -' ^ 

n=-oo m^n 1 
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CO 

im'ki 



-i /-oo 1 — _ -i 

-— 2 / dk ld k 2 a(ki,t)a*(k 2 ,t) £ e-^^-W £ 

V / J —oo _ — ,/_/Ln I I 



a+i e 



n=— oo m'/O 



(27T) 



-i /-co 00 

— / dk 1 dk 2 a(k 1 ,t)a*(k 2 ,t) y e'^+^-^Jaih) 



1 

2^ 



00 

dkidk 2 a(ki,t)a* (k 2 ,t)5(k + k 2 — ki)J a (ki] 



00 



= — / dkia(ki,t)a*{ki-k,t)J a (ki), (140) 

where J a (A;) is defined by (JH|). Multiplying the final expression (|140j) by exp(ikx), and integrating 
over from —00 to 00, we get finally 

1 /'OO 1 /'OO 

— / dke ikx — dk 1 a(ki,t)a*(kx-k,t)J a (k 1 ) = 
27r J -00 27T J -00 

-l /"CO -| f'OG 

= — / dfcio(fci,t)J Q (fci)— / dke ikx a*(ki - k,t) = 

2?r V-oo 27r J -OO 

-I /'OO 1 /"OO 

= — / dfcia(fci ) i)J a (A;i)e* a: -- / dke~ ik ' x a* (k' , t) = 

J -oo 2?r ./-oo 

1 f 00 

= 5*(x,t)— / dfcia(A:i,t)J a (A;i)e ifclX = 
2vr J^oo 



= ^ t) ( 2 C(« + 1)5(0;, t) - a a ^S(x, t) + 2 E «)) , (141) 

where a a is defined by (|121|) . and we use (|16j) . This expression is just the equation (|12U|) . 
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